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ABSTRACT 

The Advanced Camera for Surveys on board the Hubble Space Telescope has been used to obtain 
deep, high-resolution images of the intermediate-age star cluster NGC 1846 in the Large Magellanic 
Cloud. We present new color-magnitude diagrams (CMDs) based on F435W, F555W, and F814W 
imaging. We test the previously observed broad main sequence turnoff region for 'contamination' by 
field stars and (evolved) binary star systems. We find that while these impact the number of objects 
in this region, none can fully account for the large color spread. Our results therefore solidify the 
recent finding that stars in the main sequence turnoff region of this cluster have a large spread in 
color which is unrelated to measurement errors or contamination by field stars, and likely due to a 
~ 300 Myr range in the ages of cluster stars. An unbiased estimate of the stellar density distribution 
across the main sequence turnoff region shows that the spread is fairly continuous rather than strongly 
bimodal, as suggested previously. We fit the CMDs with several different sets of theoretical isochrones, 
and determine systematic uncertainties for population parameters when derived using any one set of 
isochrones. We note a degeneracy between age and [a/Fe], which can be lifted by matching the shape 
(curvature) of the full red giant branch in the CMD. We find that stars in the upper part of the main 
sequence turnoff region are more centrally concentrated than those in any other region, including more 
massive red giant branch and asymptotic giant branch stars. We consider several possible formation 
scenarios which account for the unusual features observed in the CMD of NGC 1846. 
Subject headings: galaxies: star clusters — globular clusters: general — Magellanic Clouds) 



1. INTRODUCTION 

An accurate knowledge of stellar populations of "in- 
termediate" age (« 1 - 3 Gyr) is important within the 
context of several currently hot topics in astrophysics. 
Intermediate-age stars typically dominate the emission 
observed from galaxies at high redshift (e.g., van der Wei 
et al. 2006). Furthermore, star clusters in this age range 
are critical for testing predictions of the dynamical evolu- 
tion of star clusters (e.g., Goudfrooij et al. 2007), and for 
understanding the evolution of intermediate-mass stars. 
The Large Magellanic Cloud (LMC) hosts a rich system 
of intermediate-age star clusters. The first surveys dedi- 
cated to studying properties of these clusters were based 
on integrated colors (e.g., Searle et al. 1980). These stud- 
ies led to an empirical and homogeneous age scale based 
on the so-called "5 parameter" (Elson & Fall 1985, 1988; 
Girardi et al. 1995; Pessev et al. 2008), which describes 
the position of the cluster in the integrated {U — B) ver- 
sus (B—V) color-color diagram. While a number of more 
recent studies of such intermediate-age clusters have de- 
termined ages more directly by fitting the location of 
the main sequence turnoff (MSTO) region with model 
isochrones by means of color-magnitude diagrams (CMD; 
e.g., Mighell et al. 1998; BerteUi et al. 2003; Kerber et 
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al. 2007; Mucciarelh et al. 2007; Mackey et al. 2008), 
such age determinations are still rather sparse and often 
dependent on the stellar model being used. It is impor- 
tant to obtain more CMD-based ages and metallicities 
of intermediate-age clusters and to study any systematic 
uncertainties related to the choice of any particular stel- 
lar model. 

Mackey et al. (2008) recently used deep, HST/ACS 
images to study the intermediate age LMC star cluster 
NGC 1846 (plus two other clusters from our program 
GO-10595; PI: Goudfrooij). They found evidence for a 
wide range in turnoff colors in a fairly bimodal distribu- 
tion, and concluded that these are due to two bursts of 
star formation separated by w 0.3 x 10^ yr. Recently, 
a number of star clusters, primarily massive old globu- 
lar clusters in the Milky Way have also been shown to 
have unusual color- magnitude diagrams (CMDs), sug- 
gestive of multiple stellar populations with variations 
in either age or abundance. This includes wCen (e.g., 
Norris et al. 1996; Hilker & Richtler 2000; Bedin et al. 
2004; ViUanova et al. 2007), NGC 2808 (Piotto et al. 
2007), NGC 1851 (Milone et al. 2008a), and NGC 6388 
(Piotto 2008). These globular clusters are among the 
most massive known in our Galaxy, and several are be- 
lieved to be the remnant nuclei of stripped dwarf galaxies. 
NGC 1846 is roughly an order of magnitude less massive 
than known Galactic clusters with multiple populations, 
and it seems unlikely to be the stripped nucleus of a can- 
nibalized dwarf galaxy since it is located in the outskirts 
of the LMC, which is a dwarf irregular galaxy itself. 

In this paper we conduct a more detailed investigation 
of NGC 1846 in the LMC than that presented in Mackey 
et al. (2008), including: The effects of photometric com- 
pleteness and binary star evolution, new techniques for 
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assessing background contamination, and the radial dis- 
tribution of cluster stars at different evolutionary phases. 
We also investigate how well different sets of isochrones 
and variations in [a/Fe] fare when compared with the ob- 
servations. We accomplish this by using the ePSF pho- 
tometric technique (Anderson & King 2006; Anderson et 
al. 2008a, b) to construct the cluster CMD, while Mackey 
et al. used DOLPHOT (Dolphin 2000) for their analysis. 
Our results support the general conclusion of Mackey et 
al. (2008) that there is a large spread in the main se- 
quence turnoff color for NGC 1846 which must be due 
to a range in ages, although we do find differences and 
additional relevant properties. 

The remainder of this paper is organized as follows. § 2 
presents the observations, § 3 discusses details of the stel- 
lar photometry, completeness corrections, the evaluation 
of contamination by field stars, and radial distributions 
of stars in different parts in the color-magnitude diagram 
of NGC 1846. § 4 presents our isochrone fitting analysis 
which involves different sets of theoretical isochrones, § 
5 discusses the results, and § 6 presents our main conclu- 
sions. 

2. OBSERVATIONS 

NGC 1846 was observed with HST on Jan 12th, 2006, 
using the wide-field channel (WFC) of ACS, as part of 
HST program 10595 (PI: Goudfrooij). We centered the 
cluster on one of the two CCD chips of the ACS/ WFC 
camera, so that the observations cover enough radial ex- 
tent to study variations with cluster radius. Three ex- 
posures were taken in each of the F435W, F555W, and 
F814W filters: Two long exposures of 340 s each and 
one shorter exposure to avoid saturation of the brightest 
stars (90 s, 40 s, and 8 s in F435W, F555W, and F814W 
respectively) . The two long exposures in each filter were 
spatially offset from each other by 3'.'011 in a direction 
-|-85?28 with respect to the positive X axis of the CCD 
array. This was done to move across the gap between 
the two ACS/WFC CCD chips, as well as to simplify the 
identification and removal of hot pixels. 

We used two types of products provided by the HST 
calibration pipeline, namely the fit and drz files. The 
fit files have been bias corrected, dark subtracted, and 
flatfielded. These images are still in the WFC CCD co- 
ordinate frame, which suffers from strong geometric dis- 
tortion. The drz files are composite images of exposures 
taken using the same filter, and are useful in the con- 
text of our analysis for two main reasons: (i) The drz 
files were produced using Multidrizzle (Koekemoer et al. 
2003) which corrects the images for geometric distortion 
and ties the output image to the astrometric reference 
frame of the HST guide star catalog; (ii) the photomet- 
ric calibration of the ACS cameras was performed on drz 
images (Sirianni et al. 2005, hereafter S05). On the other 
hand, the resampling of the images by the drizzle algo- 
rithm (Fruchter & Hook 2002) within the Multidrizzle 
package makes them less suited for accurate point-spread 
function (PSF)-fitting photometry. Hence, we use the 
drz files for photometric and astrometric calibration and 
the fit files for measurements. 

3. ANALYSIS 
3.1. Photometry 



Stellar photometry was performed using PSF fitting, 
using the spatially variable "effective PSF" (ePSF) 
method described in Anderson & King (2000), and tai- 
lored for ACS/WFC data by Anderson & King (2006). 
A detailed description of the application of the ePSF 
method to ACS/WFC data is given in Anderson et al. 
(2008a, b). Small temporal variations of the PSF from 
one image to another (likely due to thermal breathing 
and/or small changes in the focus) were dealt with by 
constructing a spatially constant "perturbation" PSF for 
each image which was added to the library ePSF for each 
filter. 

After fitting positions and fiuxes of all stars in all indi- 
vidual fit images using the ePiSF-fitting algorithm, they 
were transformed into the geometrically corrected (drz) 
image reference frame using the prescriptions of Ander- 
son & King (2006). We selected all stars with the ePSF 
parameter "PSF fit quality" q < 0.5 and "isolation in- 
dex" of 5. The latter parameter selects stars that have no 
brighter neighbors within a radius of 5 pixels. To further 
weed out hot pixels, cosmic rays and spurious detections 
along diffraction spikes, the geometrically corrected po- 
sitions among the three images per filter were compared. 
We selected objects with coordinates matching within a 
tolerance of 0.2 pixels in either axis, which eliminated the 
hot pixels and cosmic ray hits effectively. The photom- 
etry, at this stage, is given as instrumental magnitudes 
(—2.5 logj^Q (flux [e^])) for each individual exposure. 

Corrections for imperfect charge transfer efficiency 
(CTE) of the ACS/WFC CCDs were made following 
Kozhurina-Platais et al. (2007), using a functional form 
of CTE loss that was derived from the data itself by com- 
paring the photometry of stars from the short versus the 
long exposures. Using this tcc;linique, we find a differ- 
ent functional form for CTE loss than that published by 
Riess & Mack (2004). The latter was found to overcorrect 
the CTE loss. This difference is likely due to two main 
reasons: (i) The Riess & Mack (2004) correction formu- 
lae were derived from sparsely populated fields, whereas 
we are working with a much more crowded field. Charge 
traps in crowded fields are expected to be filled by charge 
associated with stars located at rows close to the read- 
out amplifier and hence the overall effective CTE loss 
is expected to be lower than in sparse fields (see, e.g., 
Goudfrooij et al. 2006); (ii) Riess & Mack (2004) per- 
formed aperture photometry on drz images, whereas we 
are using PSF-fitting photometry on fit files. The lat- 
ter method assigns higher weights to the central pixels of 
stars relative to those located in the wings of the PSF, 
whereas aperture photometry assigns equal weights to 
every pixel within the aperture. The relevance of the 
CTE correction is illustrated in Figs. 1 and 2. The accu- 
racy of our CTE corrections is such that the rms scatter 
of the magnitude residuals is 0.02 mag at an instrumen- 
tal magnitude of —2 (corresponding to mF435w = 23.8, 
?^F555w = 24.4, and mF8i4w = 23.5 in magnitude units 
relative to Vega). This uncertainty is smaller than the 
photometric measurement errors at those magnitudes. 
After CTE correction, a final magnitude for each object 
was determined from a weighted average from all three 
available photometric measurements in a given filter. A 
weight of was used after rejecting saturated sources 
in the long exposures. 



The Wide MSTO of NGC 1846 in the LMC 



3 



Q 2 - ITFCB 




-a -6 -4 -2 -fl -B -4 -2 

TriHtrnmpntRl Vno TTiHti-nTTi#»Tit b1 Uret 

Fig. 1. — The difference in normalized instrumental F435W mag- 
nitude (—2.5 log (e~ s~^)) derived from short vs. long exposures 
as a function of normalized instrumental F435W magnitude of the 
long exposures. The left panels show the observed values (affected 
by CTE loss), whereas the right panels show the result of the CTE 
correction. The top and bottom panels refer to stars on the WFC2 
and WFCl chips, respectively. The overplotted solid (red online) 
line illustrates the average trend in magnitude residuals. 
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Fig. 2. — The difference in normalized instrumental F435W 
magnitude derived from short vs. long exposures as a function of 
Y coordinate in the drz reference frame. The left panels show the 
observed values (affected by CTE loss), whereas the right panels 
show the result of the CTE correction. The top row of panels refer 
to stars with normalized instrumental magnitudes — 10 < mag < 
—6, the middle row does so for —6 < mag < —4, and the bottom 
row does so for —4 < mag < —2. The overplotted solid (red online) 
line illustrates the average trend in magnitude residuals for each 
chip (WFC2 for Y < 2048 and WFCl for Y > 2048). 



3.2. Photometric Calibration 

We determined final magnitudes for each star on the 
Vegamag photometric system following the procedure 
outlined in Bedin et al. (2005). Briefly, instrumental 
magnitudes measured with the ePSF method were com- 
pared with aperture photometry with a 0.5" radius of the 
same stars. Our final, calibrated magnitudes were then 
determined using: 

mfi,„.--2.5 1og,„(M£n)+ZPfi,er 



where the subscript "filter" refers to the specific filter 
that was used. The first term on the right hand side is the 
normalized instrumental magnitude measured from the 
fit images using ePSF, the second term is the zeropoint 
from 805*^, the third term is the correction from 0'.'5 to 
an infinite aperture (taken from table 1 in Bedin et al. 
2005), and the fourth term is the empirically determined 
difference between our ePSF-fitting photometry and the 
aperture photometry performed for bright isolated stars 
in the image. 

3.3. Color-Magnitude Diagrams 

Color-magnitude diagrams (CMDs) for NGC 1846 
were created for both mF435w vs. tof435W — ™f814W and 
'T^Fsssw vs. mF555w ~ '71F814W- AH stars selected as de- 
scribed in § 3.1 above are plotted in panels (a) and (b) of 
Fig. 3. To decrease contamination by field stars, we also 
plot CMDs including only stars within 25 arcsec from 
the center of NGC 1846 (which has reference pixel coor- 
dinates (xc, yc) = (2225, 3096), see § 3.5 below) in panels 
(c) and (d). There are several features clearly visible in 
the CMD of NGC 1846 shown in Fig. 3: 

1 . The main sequence is well defined and extends over 
about 6 mag in tof435W and tofsssw, and is ac- 
companied by an obvious sequence of unresolved 
binary stars (running along the main sequence but 
slightly brightward and redward). The main se- 
quence shows a slight change in slope near tof435w 
K, 22.5, flagging the transition between radiative 
and convective stellar cores. Fig. 3 shows that the 
MSTO region is fairly wide, whereas the MS of sin- 
gle stars fainter than the turnoff is well-defined and 
much narrower (see especially panel (e)). This con- 
stitutes strong evidence for the presence of more 
than one population, which is discussed in more 
detail below. 

2. The subgiant branch (SGB) is quite well defined 
and seems to emerge mainly from the upper end of 
the MSTO region. 

3. The red giant branch (RGB) is well defined and 
quite narrow (see especially panels (a) and (e) of 
Fig. 3). This strongly suggests that there is little 
differential reddening across NGC 1846. Note that 
the significant width in color of the MSTO region is 
not reflected in the RGB. This argues that the pop- 
ulations that make up the MSTO share the same 
metallicity (in terms of [Fe/H]). 

The data allow one to clearly identify the so-called 
"RGB bump" at mF435w = 20.0, mF435w— "^F814W 
= 2.1, mpsssw = 19.0, and mpsssw - "^fbmw = 
1.1. This feature was first predicted by Thomas 
(1967) and Iben (1968), and corresponds to the 
period during Hydrogen shell burning when the 
outward-moving shell encounters extra Hydrogen 
fuel left behind during the period when the star 
had a convective envelope. The "bump" occurs as 
the star adjusts its structure as it consumes this 
extra fuel, briefly reversing its path through the 

^ The SOS zeropoints were recently updated 
for different observing dates, and are available at 
http: / / www. stsci. edu /hst /acs /analysis / zeropoints . 
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Fig. 3. — CMDs for the HST/ACS images of NGC 1846. The top panels (a) and (b) show mF435w vs. mF43r)W — "IF814W and "IF555W 
vs. mF555W — "'FS14W for all 23,652 detected sources that passed our selection criteria. The middle panels (c) and (d) do so for a region 
at radii 125" -190" from the center of NGC 1846 (near the corners of the HST/ACS image), and the bottom panels (e) and (f) do so for 
the selected sources within 25" from the center of NGC 1846, a region with the same area as that shown in panels (c) and (d). 



CMD before resuming its ascent of the RGB (see 
also, e.g., Fusi Pecci et al. 1990; Denissenkov & 
VandenBerg 2003). 

4. The red clump (RC), sometimes called the Helium 
clump (since it represents core Helium burning), is 
located at tof435w — 20.1, mF435w ~ mpsiAw — 
1.9, mF555W - 19.2, and mF555W - mpsww - 1-0. 

5. The asymptotic giant branch (AGB) clump can be 
seen at mF435w — 19.5, mF435w — "^f814W — 2.2, 
TOF555W - 18.4, andmF555W-"^F8i4w - 1-2. This 
feature corresponds to the onset of Helium shell 
burning. 

The wide MSTO region in NGC 1846 was previously 
reported by Mackey & Broby Nielsen (2007) and Mackey 
et al. (2008), who interpreted this feature as due to the 
presence of two stellar populations of identical metal- 
licity but different age. Our photometry shows some 
differences with that of Mackey and coworkers that po- 



tentially affects the interpretation of how this cluster 
formed. These differences are discussed further in § 5. 

3.4. Completeness Corrections 

Density distributions of stars in crowded images suffer 
from incompleteness. We use the standard technique of 
adding artificial stars to the images and running them 
through the photometric pipeline in order to quantify 
incompleteness. Specifically we use a program written 
by Jay Anderson (see Anderson et al. 2008a) that uses 
the best-fit arrays of ePSFs for each individual image, 
repeatedly adding artificial stars covering the magnitude 
and color ranges found in the CMDs. The overall radial 
distribution of the inserted artificial stars followed that 
of the stars in the image. Only 5% of the total number of 
stars selected in § 3.1 above were added per run, so as not 
to increase the degree of crowding in the images signifi- 
cantly. After inserting the artificial stars, the ePSF mea- 
surement and selection procedures were applied again to 
the image. An inserted star was considered recovered 
if (i) the input and output positions agreed to within 
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0.3 pixel, and (ii) the input and output fluxes agreed to 
within 0.5 mag in every filter. Completeness fractions of 
stars as a function of magnitude and distance from the 
cluster center are shown in Fig. 5 for stars in the cluster 
sequences, as depicted in Fig. 4. 

Completeness fractions were assigned to every individ- 
ual star in a given CMD by fitting the completeness frac- 
tions as a function /comp('7i, r) of magnitude m and ra- 
dius r: 



ellipse and e its ellipticity. The radial surface number 
density profile was fit with a King (1962) model com- 
bined with a constant background level: 



1 



(\\r) (III - ;/)um(/-)) 



(1) 



where a{r) is a free scaling parameter and mumir) is the 
magnitude where completeness is 50%. The functional 
forms of both a{r) and miiin(r) were chosen as third- 
order polynomials, fit to the data using 10 bins in radius. 
The overall rms of the fit to the completeness data was 
0.03 for both the mF4:55w vs. mF435w — tof814W and the 

TOF555W vs. rnF555W - TOF814W CMDs. 

3.5. Radial Distributions 

The wide main sequence turn-off region in NGC 1846 
suggests the existence of more than one simple popu- 
lation. If so, it is important to find out whether the 
different populations have intrinsically different spatial 
distributions, i.e., differences beyond those that can be 
expected for simple (coeval) stellar populations due to 
dynamical evolution. A well-known example is mass 
segregation due to dynamical friction which slows down 
stars on a time scale that is inversely proportional to the 
mass of the star (e.g., Saslaw 1985; Spitzer 1987; Mey- 
lan & Heggie 1997) so that massive stars have a more 
centrally concentrated distribution over time than less 
massive ones. With this in mind, we determined regions 
in the CMD that are strongly dominated by stars be- 
longing to the cluster rather than to the underlying field 
in the LMC. 

3.5.1. Surface Density Distribution 

We first determine the projected number density of 
stars, taking into account the completeness corrections. 
The cluster center was determined to be at reference co- 
ordinate {xc,yc) = (2225,3096) with an uncertainty of 
± 5 pixels in either axis. This was done by creating a 
2-D histogram of the pixel coordinates of stars using a 
bin size of 50 x 50 pixels (i.e., 2'.'5 x 2'.'5), and then cal- 
culating the center using a 2-D gaussian fit to an image 
constructed from the surface number density values in 
the 2-D histogram. Note that this method avoids biases 
related to the presence of bright stars near the center. 
The ellipticity of NGC 1846 was found to be 0.12 ± 0.02, 
which was derived by running the task ellipse within 
IRAF/stsdas^ on the surface number density image men- 
tioned above. The area sampled by the ACS image was 
divided in 12 concentric elliptical annuli, centered on the 
cluster center. For annuli with radii larger than ^ 850 
pixels (= 42'.'5), care was taken to account for the lim- 
ited azimuthal coverage of the cluster in the image. Radii 
are expressed in terms of the "equivalent" radius of the 
ellipse, r — a yT^e where a is the semimajor axis of the 

^ STSDAS is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA 



n(r) = no 
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+ bkg (2) 



where no is the central surface number density, Tc is the 
core radius, and c = (rt/rc) is the concentration index 
{rt being the tidal radius). The best-fit King model was 
selected using a minimization routine, and is shown 
in Fig. 6 as well as the individual surface number density 
values. 

3.5.2. Selection of Cluster-Dominated Regions on the 

CMD 

The next step in establishing regions on the CMD that 
are strongly dominated by cluster stars was done sta- 
tistically, by comparing star surface number densities in 
selection boxes on the CMD from two radial ranges: "In- 
ner" stars with log(r) < 1.5 versus "outer" stars with 
2.0 < log(r) < 2.2 (with radius r in arcsec). Selection 
boxes for which the "inner/outer" surface number den- 
sity ratio (after completeness correction) exceeded the 
value given by the best-fit King model to all stars in the 
ACS image were tagged as dominated by cluster stars. 
The selection boxes themselves were created for each evo- 
lutionary sequence in the CMD (MS, RGB, He clump, 
AGB), plus (bigger) boxes in the less populated areas 
of the CMD. For each evolutionary sequence, a ridgeline 
was established by means of a spline fit to carefully se- 
lected reference points along the sequence in question. 
After subtracting the color of the ridgeline from those of 
the stars in that sequence, selection boxes were divided 
into bins of magnitude and color. The bin size for the 
MS region was chosen adaptively depending on the local 
number density in the CMD, with a minimum of 0.2 mag 
in magnitude and color. For the RGB and AGB, the bin 
size in color was chosen to properly separate the two se- 
quences. The selection boxes that did not correspond to 
the evolutionary sequences are depicted as such on Fig. 
7 with dashed lines. 

After discarding the regions not dominated by cluster 
stars, we derived radial distributions of the following re- 
gions on the CMD, all of which were found to have less 
than 10% contamination by field stars: (i) the lower MS, 
meaning the MS below the turn-off of the field stars; (ii) 
the binary sequence associated with the lower MS; (Hi) 
the upper MS, meaning the MS below the turn-off area 
but above the turn-off of the field stars; (iv) the binary 
sequence associated with the upper MS; (v) the 'lower 
half of the MSTO region; (vi) the 'upper half of the 
MSTO region; (vii) the RGB and AGB stars. Fig. 7 de- 
picts these regions on the CMD with solid lines, while 
the completeness-corrected radial distributions of stars 
in these regions are shown in Fig. 8. Only stars with 
completeness fractions > 50% were used for this exercise. 

Fig. 8 shows several items of interest. First of all, a 
comparison of the open squares with the filled squares 
in the left panel of Fig. 8 reveals that the lower MS is 
more heavily contaminated by field stars than the upper 
MS. This could have been expected given the presence 
of a wide turn-off near tof435w = 23.0 and mpsssw = 
22.5 (cf. Figs. 3a, b), which is clearly due to field stars. A 
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Fig. 4. — CMDs of stars in regions in color-magnitude space that have less than 10% of field star contamination (see text in § 3.3). 
Left panel: mF435W vs. mF435w — »tif814W- Right panel: mpsssw vs. mpsssw — "1P814W- The open red circles in the CMDs depict the 
magnitudes and colors used for the completeness fractions shown in Fig. 5 (see § 3.4). 
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Fig. 5. — The completeness fraction as function of magnitude and radius. Open squares and solid lines represent data at a radius of 5", 
open circles and dashed lines represent data at a radius of 35", and open triangles and dotted lines represent data at a radius of 125". 
Error bars depict the standard deviations among 50 independent runs of the artificial star tests. See description in § 3.4. 
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Fig. 6. — Radial surface density profile of NGC 1846. The points 
represent observed values. The dashed line represents the best-fit 
King model (cf. Eq. 2) whose parameters are shown in the legend. 



similar effect is seen for the binary MS sequence (com- 
pare open circles with filled circles in the right panel of 
Fig. 8). Second, a comparison of the radial distribution 
of the upper MS stars with that of the bright RGB/AGB 



stars (filled squares and filled triangles in the left panel of 
Fig. 8, respectively) shows that mass segregation is not 
particularly strong in NGC 1846. There is a hint of the 
effect in the inner ~ 20" (equivalent to ~5 pc), where 
the surface number density of RGB/AGB stars contin- 
ues to rise towards the center at a slightly stronger rate 
than that of the upper MS stars. Last but not least, 
the right panel of Fig. 8 shows that the stars in the up- 
per (brighter) half of the MSTO region are significantly 
more centrally concentrated in NGC 1846 than those in 
the lower (fainter) half of the MSTO. In fact, the stars in 
the upper half of the MSTO show a central concentration 
that is similar to (or even stronger than) that of stars in 
the upper RGB/AGB. The very small difference in the 
luminosity and mass (i.e., < 5%) of stars in the upper vs. 
the lower halves of the MSTO region strongly suggests 
that biases and mass segregation are not the cause of this 
difference. We conclude that the upper and lower MSTO 
regions correspond to intrinsically different populations 
which may well have undergone different amounts of vio- 
lent relaxation during their collapse and/or different dy- 
namical evolution effects. We discuss potential origins 
for this finding in § 5. 

4. ISOCHRONE FITTING 
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Fig. 7. — A copy of Fig. 3a showing regions of interest. The 
red solid lines encompass regions in color-magnitude space that 
have less than 10% of field star contamination and for which radial 
distributions are shown and described in Fig. 8. The blue dashed 
lines encompass regions found to be heavily contaminated by field 
stars. See discussion in § 3.5.2. 



We fit isochrones to the CMDs of NGC 1846 to deter- 
mine the age, metalhcity, [a/Fe] ratios, and their associ- 
ated uncertainties. Three sets of stehar models with pre- 
dictions computed for the ACS/WFC filter system were 
used: Padova isochrones (Marigo et al. 2008; Girardi et 
al. 2008), Teramo isochrones (sometimes referred to as 
BaSTI isochrones; Pietrinferni et al. 2004, 2006), and 
Dartmouth isochrones (Dotter et al. 2008). 

Padova isochrones: We use the default models which 
involve scaled solar abundance ratios (i.e., [a/Fe] = 0.0) 
and which include some degree of convective overshoot- 
ing (see Girardi et al. 2000). Using the web interface of 
the Padova team'', we construct a grid of isochrones that 
covers the ages 0.3 < T[Gyr] < 3.0 (where r is the age) 
with a step of At = 0.1 Gyr and metallicities Z = 0.001, 
0.002, 0.004, 0.006, 0.008, 0.01, 0.02, and 0.03. 

Teramo isochrones: We use that team's web site^ 
to construct grids of isochrones that cover the same ages 
and metallicities as for the Padova models mentioned 
above. We use the Teramo isochrones with [a/Fe] = 
0.0. To allow an assessment of the influence of different 
amounts of convective overshooting, we construct grids 
for both the "canonical" models which do not include 
overshooting and the "non-canonical" models which do. 

Dartmouth isochrones: We use the full grid avail- 
able from their web site^ which covers the ages 0.25 < 
r[Gyr] < 1.0 with At = 0.05 Gyr and 1.0 < r[Gyr] < 5.0 
with At = 0.25 Gyr, metallicities [Fe/H] = -2.5, -2.0, 
-1.5, -1.0, -0.5, 0.0, +0.3, and +0.5, and [a/Fe] = 
-0.2, 0.0, +0.2, +0.4, +0.6, and +0.8. 

4.1. Fitting Method 

http: //stev . oapd. inaf . it/cgi-bln/cmd 
* http : //albione . oa-teramo . inaf . it 

® http: //stellar .dartmouth . edu/~models/complete . html 



Isochrone fitting is performed as follows. We first es- 
tablish parameters that involve pairs of fiducial points on 
the CMD that are: (i) relatively easy to measure or de- 
termine from both the data and the isochrone tables, (ii) 
sensitive to at least one population parameter such as age 
or metalhcity, and (Hi) independent on the distance and 
foreground reddening of the cluster. This allows us to 
reduce the large number of free parameters, which is al- 
ways a concern when fitting isochrones to cluster CMDs. 
After extensive experimentation, we select the following 
parameters: 

1. The difference in magnitude between the MSTO 
and the RGB bump, defined as AB^Jg^ and 
^^rgbb' B and V filters, respectively. The 
MSTO is defined as the point where a polynomial 
fit to the stars (or the isochrone entries) near the 
turn-off is vertical in the CMD. For the location of 
the RGB bump on the CMD we simply calculate 
the mean magnitude and color of stars in a box cen- 
tered on the RGB bump by eye. For the isochrones, 
we define the location of the RGB bump as the 
average of the magnitudes and colors of isochrone 
RGB entries between the two masses at which the 
magnitudes and colors "turn around" in direction 
on the CMD with increasing stellar mass. 

2. The difference in color between the MSTO and the 
RGB bump, defined as A{B - /)^Jgg and A{V - 

nMSTO 
-'7RGBB- 

3. The slope of the RGB. This was evaluated using the 
(mean) color of the RGB stars at two fiducial mag- 
nitudes, namely at mRGBB + 1 and mRQBB — 0.75. 
The former magnitude was chosen to represent a 
point intermediate between the RGB bump and the 
lower end of the RGB; the latter magnitude was 
chosen to avoid issues related to confusing RGB 
with AGB stars on the CMD. The mean colors 
of the RGB stars were derived from the CMD by 
means of a polynomial fit to the RGB star positions 
in the CMD. In case of the isochrones, the colors 
were derived by means of linear interpolation be- 
tween isochrone table entries. 

The sensitivity of these parameters to population pa- 
rameters is illustrated in Figs. 9-11 for the age range 1 - 
3 Gyr. These plots show the measurements of ABp^Qgg , 

A{B — /)rgbB' ^^'^ ^^'^ RGB slope as functions of age 
and metalhcity (and [a/Fe] in case of the Dartmouth 
isochrones) for isochrones in a given family, compared 
to the values measured in the CMD of NGC 1846. As 
might have been expected, the RGB slope is highly sen- 
sitive to metalhcity and almost independent of age in 
the age range studied here. Figs. 9-11 also show that 
ABpQBB is dependent on both age and [Fe/H], and that 
A{B — /)rgbb primarily dependent on age, although 
there are also dependences on metalhcity and [a/Fe] (see 
Fig. 11 for the latter) among metallicities [Fe/H] < —1.0 
(the exact value of [Fe/H] depending on the isochrone 
family used). 

We then select all isochrones (within each family) for 
which the values of these three parameters lie within 2 
a of the measurement uncertainty of those parameters 
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Fig. 8. — Radial surface number density distributions of stars in different regions of the CMD, as depicted in Figure 7 with red solid 
lines. Left panel: Solid circles: all stars in the CMD. Open circles: stars in the regions of the CMD determined to have less than 10% 
contamination by field stars. Open squares: lower MS stars. Filled squares: upper MS stars. Filled triangles: RGB and AGB stars. Right 
panel: Solid circles: stars in lower binary sequence. Open circles: stars in upper binary sequence. Filled squares: lower MSTO region. Filled 
triangles: upper MSTO region. The absolute surface number density values on the Y axis refer to the solid circles on the left panel ("all 
stars"). The surface number densities of the other star types are normalized to that of "all stars" at the outermost radial bin. Error bars 
are only drawn if they are larger than the symbol size. See discussion in § 3.5.2. 
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Fig. 9. — Dependence of three distance— and reddening-independent parameters involving fiducial points on the mF435w vs. mp435w — 
'nF8i4W CMD on age and metallicity for the Padova isochrones (Marigo et al. 2008). Top panels: Slope of the RGB as defined in § 4.1. 
Middle panels: Brgbb ~ ^MSTO- Bottom panels: (B — /)rgbb ~ (^ ~ ^)msto ■ Error bars in the left and right panels reflect the variation 
of the parameter values among the isochrones with different metallicities and ages, respectively. The dashed lines in each panel represent 
the measurements of these parameters on the CMD. The upper and lower dashed lines in the middle and bottom panels refer to the upper 
and lower halves of the MSTO region, respectively. 



on the CMDs. This yielded 6-15 isochrones depend- 
ing on the isochrone family. For these isochrones, we 
then find the best-fit values for distance (m — M)q and 
foreground reddening Ay by means of a least squares 
fitting program that runs through a grid of (m — M)q 
and Ay values to compare the magnitudes and colors of 
the MSTO and the RGB bump of NGC 1846 with those 



measured from the isochrones. The grid encompasses 
18.0 < (m - M)a < 19.0 with a step of 0.01 mag and 
0.0 < Ay < 0.5 with a step of 0.01 mag. For the redden- 
ing we use Af435w = 1.351 Ay, Ap^^^yj = 1.026 Ay, and 
Af814W — 0.586 Ay. These values were derived using the 
ACS /WFC filter curves available in the synphot package 
of STSDAS and the reddening law of CardeUi et al. (1989) 
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assuming Ry — 3.1. The resulting values of (m — M)o 
and Av derived from the positions of the MSTO and 
RGB bump are averaged using inverse variance weights. 

Finally, the isochrones were overplotted onto the 
CMDs for visual examination. This revealed that there 
was a small but systematic offset in [Fe/H] between the 
best-fit isochrones for tof435W vs. mF435w — '71f814w and 
"T-Fsssw vs. mF555w — WF814W in the sense that [Fe/H] 
was always higher for the isochrone fit to mpsssw vs. 



'7^F555W — "T-F814W than to mF435W vs. TOF435W — "^F814W- 

This effect was present for all isochrone families used, 
although it was least significant for the Dartmouth 
isochrones where the same [Fe/H] yielded the best fit 
to both CMDs; in the latter case the visual check merely 
showed a slight mismatch of the RGB on the mF555W 
vs. TOF555W — ™f814w CMD (see Fig. 15). We suggest 
that it is likely that this effect is related to the fact that 
Kurucz model spectra of cool RGB stars contain more 
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'"F435W " '"F814W '"F555W " '"F8UW 



Fig. 12. — Best-fit isochrones from tfie Padova family (cf. Table 1) superposed onto the CMDs shown in Fig. 4. The solid line is the best 
fit to the upper half of the MSTO region, while the dashed line is the best fit to the lower half of the MSTO region. 
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Fig. 13. — Same as Fig. 12, but now for the Teramo isochrones that include convective overshooting. 



flux in the range ~ 5000 - 6500 A (including much of the 
V band) than empirical star spectra from the Pickles 
(1998) library at the same stellar type, as recently re- 
ported by Maraston et al. (2008, and private communi- 
cation). Model SEDs therefore have bluer V — I colors 
than observed for RGB stars, consistent with what we 
see (see Figs. 12-15). While a detailed investigation of 
the potential cause(s) of this effect is beyond the scope of 
this paper, we note that the Dartmouth isochrones used 
the Phoenix model atmospheres (e.g., Hauschildt et al. 
1999) to derive their Tes - color relations, whereas the 
Padova and Teramo isochrones used the ATLAS9 mod- 
els of R. L. Kurucz (e.g., CasteUi & Kurucz 2003). The 
use of different model atmospheres can at least partly 
explain this effect since the Phoenix models include 550 
million molecular transitions, while the ATLAS9 mod- 
els only include a few molecular lines in a semi-empirical 
manner. For example, the molecular band of MgH near 
5175 A which is in the V band and strong in cool stars 
is not included in the ATLAS9 models. Because of this 
effect, we focus on the tof435W vs. tof435w — wf814W 
CMD in terms of deriving population parameters and 
use the tofsssw vs. mFsssw — ™f814W CMD mainly for 



comparison purposes during the remaining analysis. 

In order to find the difference or spread in age that 
can explain the wide morphology of the MSTO region of 
NGC 1846, the fitting process mentioned above was con- 
ducted separately for the MSTO's measured from the 
upper and the lower halves of the MSTO region. Dur- 
ing this process, we noticed that the age step At = 0.25 
Gyr used for the grid of Dartmouth isochrones is too 
coarse for this purpose. Hence, we created additional 
Dartmouth isochrones covering a finer grid in age and 
metallicity around the values found by the procedure 
mentioned above, using the web interface of the Dart- 
mouth team web site. 

The best-fit isochrones and their population parame- 
ters (both for the upper and the lower half of the MSTO 
region) are shown in Figs. 12-15, superposed onto the 
CMDs. While the best-fit isochrones of each model fam- 
ily generally match the various stellar sequences well, 
there are some differences between the model families. 
One such difference among the fits is seen in the MSTO 
region, where the Teramo model without overshooting 
provides a poorer fit than the other three models (which 
all do incorporate some level of overshooting). This is 





briefly discussed in § 4.2 below. Another, more signif- 
icant difference among the isochrone fits is seen on the 
upper RGB, where the Padova and Teramo isochrones 
appear bluer than the observed stars whereas the Dart- 
mouth isochrone is an excellent fit to the full RGB. This 
difference is further discussed in § 4.3 below. 

The population parameters of the best-fit isochrones 
are listed in Table 1 for each isochrone family. Table 1 
also lists the final adopted parameters for NGC 1846, 
including their uncertainties that reflect both errors as- 
sociated with the isochrone fitting and systematic uncer- 
tainties related to the use of the different stellar models. 
The latter are discussed further in § 4.4 below. 

4.2. Influence of Overshooting Parameter 

Table 1 shows that the inclusion of convective over- 
shooting in isochrone models in the age range appropri- 
ate for GCs like NGC 1846 has a significant effect on 
the derived age and [Fe/H]. In the case of the Teramo 
isochrones, the inclusion of overshooting results in an in- 
crease in the derived age of ~ 15% and a decrease in the 
derived [Fe/H] of ~ 17%. 

Focusing on the morphology of the MSTO region in 



the CMDs, a comparison of the behavior of the best-fit 
Teramo isochrones that include overshooting with those 
that do not (see Figs. 13 and 14) suggests that over- 
shooting is occurring in NGC 1846. Hence we do not use 
the parameters of the best-fit Teramo isochrones that do 
not include overshooting in evaluating the final adopted 
values of the population parameters of NGC 1846. 

4.3. Influence of [a/Fe] Abundance Ratio 

We assess the infiuence of non-solar [a/Fe] abundances 
on the derived age and metallicity by comparing our 
CMD with Dartmouth ischrones for different values of 
[a/Fe]. The result is illustrated in Fig. 16, which shows 
the best-fit Dartmouth isochrones for [a/Fe] = —0.2, 
0.0, -1-0.2, and +0.4 superposed onto the tof435w vs. 
"^F435W— "^F8i4W CMD of NGC 1846. The fits were done 
to the upper half of the MSTO region in this case. Note 
that all four isochrones fit the MSTO and RGB bump lo- 
cations well, which is likely due to the fitting method we 
used (as described above in § 4). However, the detailed 
fits along the RGB differ significantly from one [a/Fe] 
value to another. In particular, there is a clear correla- 
tion between [a/Fe] and the slope or the curvature of the 
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TABLE 1 

Best-fit ages, metallicities, distances, reddening values, and [«/Fe] ratio for NGC 1846 using different isochrone families. 

See discussion in § 4.1 and 4.4. 



Isochrone Age Age pe/H] (m — M)o Ay rms error* 



Family (upper) (lower) (upper) (lower) 



Padova 


1.60 


± 


0.03 


1.90 


± 


0.03 


-0.40 


± 


0.05 


(0.0) 


18.42 


± 


0.02 


0.07 


± 


0.02 


0.04 


0.06 


Toramo w/ OS'' 


1.50 


± 


0.03 


1.75 


± 


0.03 


-0.40 


± 


0.05 


(0.0) 


18.50 


± 


0.02 


0.08 


± 


0.01 


0.07 


0.03 


Tcramo w/o OS^ 


1.30 


± 


0.05 


1.60 


± 


0.05 


-0.32 


± 


0.05 


(0.0) 


18.43 


± 


0.02 


0.03 


± 


0.02 


0.03 


0.02 


Dartmouth 


1.75 


± 


0.03 


2.05 


± 


0.03 


-0.50 


± 


0.05 


+0.2 ± 0.1 


18.41 


± 


0.02 


0.09 


± 


0.02 


0.02 


0.02 


Dartmouth^ 


1.95 


± 


0.03 


2.25 


± 


0.03 


-0.47 


± 


0.05 


(0.0) 


18.40 


± 


0.02 


0.11 


± 


0.02 


0.03 


0.03 



Adopted Values 1.7 ±0.2 2.0 ± 0.2 -0.5 ± 0.1 +0.2 ± 0.1 18.45 ± 0.05 0.09 ± 0.03 



*RMS error of least-square fit to reference points on the mF43Bw 
vs. mF436w ~ "IF814W CMD, see discussion in § 4.1 
''Teramo isochrones that include overshooting (OS) 
"^Teramo isochrones that do not include overshooting 
'^Dartmouth isochrones with [a/Fe] = 0.0 only 



RGB in the sense that larger [a/Fe] values give 'flatter' 
slopes and stronger curvature for the RGB. 

Furthermore, this comparison shows that larger values 
of [a/Fe] result in younger fitted ages, and hence there 
is a degeneracy between age and [a/Fe] if one doesn't 
take the detailed morphology of the RGB into account. 
The amplitude of this effect is roughly 20% in age for 
a difference in [a/Fc] of 0.2 dex (see Fig. 16). Finally, 
note that the effect of the chosen value of [a/Fe] to the 
best-fit [Fe/H] is small within the context of our fitting 
procedure: A [Fe/H] /A [a/Fe] ~ —0.15. Note however 
that the effect of non- solar [a/Fe] to determinations of 
[Fe/H] can be significantly larger for analyses that aim to 
determine [Fe/H] values for more distant stellar systems 
where only the upper end of the RGB is available for 
isochrone fitting. This effect will be quantified in detail 
in a separate paper. 

As Fig. 16 shows, the best fit to the full RGB is 
achieved using the isochrone with [a/Fe] = -1-0.2. Since 
this fit is clearly better than any isochrone that uses solar 
abundance ratios in any of the three families, we adopt 
[a/Fe] = +0.2 for NGC 1846. 

4.4. Systematic Uncertainties in Derived Properties 

We quantify systematic uncertainties in derived age, 
[Fe/H], distance, and reddening by comparing our best- 
fit results from each set of isochrones with convcctivc 
overshooting and [a/Fe] = 0.0, as compiled in Table 1. 
These results give systematic uncertainties of: ±0.24 
Gyr in age (~ 15%), ± 0.05 dex in [Fe/H] (~ 12%), ± 0.05 
mag in (m — M)o (~5% in linear distance), and ±0.03 
mag in Ay (— 30%). We suggest that these values repre- 
sent typical systematic uncertainties associated with the 
determination of population parameters of intermediate- 
age star clusters from rnF435w vs. mF435w — 'tif814W 
CMD fitting by isochrones of any given stellar model. 

Finally, we note that the properties of the best-fit 
isochrones arc consistent with an estimate of the mass 
of AGB stars in NGC 1846 as derived from their pul- 
sational properties: Lebzelter & Wood (2007) estimated 
■Magb — 1-8 Mq with an uncertainty of ~ 0.2 Mq 
(Lebzelter, private communication). The formal AGB 
star masses for the best-fit isochrones in Table 1 all fall 
between 1.70 and 1.77 Mq, consistent with the indepen- 



dent estimate of Lebzelter & Wood. 

5. DISCUSSION 
5.1. Morphology of the MSTO Region 

As already mentioned in §§ 3.3 and 3.5.2, the MSTO 
region is broader than can be explained by photomet- 
ric errors, and the upper (brighter) MSTO region has a 
significantly more centrally concentrated radial distribu- 
tion than the lower one. This finding strongly suggests 
the presence of multiple (at least two) stellar populations 
in NGC 1846. To constrain this picture further in terms 
of whether or not two distinct populations arc sufficient 
(or necessary) to produce the observed CMD, we conduct 
simulations of a synthetic cluster with the properties im- 
plied by the isochrone fitting in § 4. For this purpose we 
define the null hypothesis such that the upper and lower 
halves of the MSTO region represent two distinct simple 
stellar populations with properties given in Table 1 for 
the best-fit Dartmouth isochrones. 

We simulate a simple stellar population (SSP) with a 
given age and chemical composition by populating an 
isochrone with stars randomly drawn from a Salpeter 
mass function between 0.1 Mq and the RGB-tip mass, 
which is at 1.62 Mq for the best-fit Dartmouth isochrone 
(see 4th entry of Table 1). The total number of stars is 
normalized to the observed number of stars brighter than 
the 50% completeness limit. To a fraction of this sam- 
ple of stars we add unresolved secondary binary compo- 
nents that are drawn from the same mass function, hence 
assuming a fiat primary-to-secondary mass ratio distri- 
bution for binary stars. We note that, although we are 
primarily investigating the morphology of the MSTO re- 
gion, we fully model the stellar evolution until the onset 
of He-shell burning and also add secondary binary com- 
ponents that were evolved past the tip of the RGB. Their 
evolution is modeled using the He-core burning evolu- 
tionary tracks assuming an average mass-loss of 0.15 Mq 
along the RGB and a stellar mass dispersion on the zero- 
age horizontal branch of 0.06 Mq (sec Lec 1991; Dottcr 
2008, for details). We use the width of the upper main 
sequence, i.e. the part brighter than the the turn-off of 
the background stellar populaticm, to fix the binary star 
fraction at 15%. We estimate the internal systematic 
uncertainty in binary fraction to ±5% and defer a more 
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detailed discussion of the parameter degeneracies, in par- 
ticular binarity vs. mass fractions of stellar generations 
(i.e., the star formation history of the cluster), to a future 
paper. For the purposes of this work the results don't 
change significantly within ~ 10% of the binary fraction. 
Finally, we add photometric errors that are modeled us- 
ing the distribution of photometric uncertainties of our 
observations to the sample of artificial stars. 

We compare the observed MSTO morphology of 
NGC 1846 with those of the simulated CMDs by con- 
structing a parallelogram in the CMD with the following 
characteristics: (i) One axis is approximately parallel to 
both isochrones, (ii) the other axis is approximately per- 
pendicular to both isochrones, and (in) it is located in 
a region of the JVISTO where the split between the two 
isochrones is relatively evident. The comparison between 
data and simulations is then done by transforming the 
(mF435w — "^F8i4w, wf435w) Coordinates of the stars in 
the CMD into the reference coordinate frame defined by 
the two axes of the parallelogram, and then considering 
the distributions of reference coordinates of the stars in 
the direction perpendicular to the isochrones. 

The procedure mentioned above is illustrated in 
Fig. 17. The top panels (a, c, and e) show the CMDs 
of the data and two simulations with relevant fractions 
of stars in the two distinct SSPs (in black dots) along 
with the isochrones (in orange and green) and the par- 
allelogram mentioned above (in blue, with the reference 
axis along the isochrones in red). The bottom panels 
(b, d, and f) show the distributions of the stars' coor- 
dinates in the direction perpendicular to the isochrones 
for the CMDs shown in the corresponding top panels. 
To avoid issues related to the use of finite histogram bin 
sizes in the presence of relatively few data points, the 
coordinate distributions were calculated by means of the 
non-parametric Epanechnikov-kernel probability density 
function (Silverman 1986) for all objects in the parallel- 
ograms. In the case of the observed CMD (i.e., panels a 
and b), this was done both for stars within a radius of 



42'.'5 from the cluster center and for the "background re- 
gion" for which the CMD was shown in Fig. 3e. The area 
of the probability density function of the background re- 
gion was then scaled to that of the inner 42'.'5 radius by 
area on the ACS image, and the true probability density 
function of the inner 42'.' 5 radius was derived by statisti- 
cal subtraction of the background region (see panel b). 

As panel (b) of Fig. 17 shows, the distribution of ob- 
served stars in the parallelogram peaks at the location of 
the "younger" isochrone and then declines more or less 
uniformly towards higher coordinate values, i.e., older 
ages. A comparison with panels (d) and (f) shows that 
the distribution of the observed stars is clearly more uni- 
form than those of the two-SSP simulations. Two-sample 
Kolmogorov-Smirnov (K-S) tests confirm this visual im- 
pression: The p values of the K-S tests to compare the 
distribution of the simulated CMDs of two SSPs with 
various mass fractions of the younger population against 
the data never exceed 11% (see Table 2). It thus seems 
fair to conclude from these data that the observed CMD 
of NGC 1846 can be better explained by a population 
with a distribution of ages rather than by two discrete 
SSPs as suggested by Mackey et al. (2008) and Milone 
et al. (2008b). 

5.2. On The Formation of NGC 1846 

The results described above naturally raise the ques- 
tion: How did an intermediate-age star cluster like 
NGC 1846 with its wide MSTO form? Is its formation re- 
lated in some way to the ancient globular clusters in our 
own Galaxy for which evidence suggests the presence of 
multiple stellar populations? Based on the results of the 
experiments described above, we assume that the wide 
MSTO in NGC 1846 is due to a spread in the ages of 
stars found within this cluster. There are two general 
situations which can lead to the presence of stars with a 
wide range of ages within a single cluster: one involving 
stars which formed originally outside of the cluster and 
then were accreted in some way, and another where stars 
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the best-fit Dartmouth isochrones to the upper and lower halves of the MSTO region, as shown in Fig. 15. Panel (b): Distribution of stars 
in the parallelogram shown in the top panel, derived using a non-parametric density estimator (see text for details). The black dashed 
line represents all stars, the blue dashed line represents stars in the background area (scaled to the full area of the ACS image), and the 
solid line represents "all stars minus background". The red arrow in panel (a) indicates the positive direction of the X axis of panel (b), 
while the red line indicates its zero point. Panels (c) and (d): Same as panels (a) and (b), respectively, but for the simulation of two SSPs 
including binary stars as described in § 5.1, with 50% of the stars belonging to the 1.75 Gyr old SSP. Panels (e) and (f): Same as panels 
(c) and (d), respectively, but now for the case where 60% of the stars belong to the 1.75 Gyr old SSP. 



TABLE 2 

Results of two-sample Kolmogorov-Smirnov tests 
regarding the morphology of the msto region. null 

hypothesis is that THE DATA OF NGC 1846 AND THE TWO-SSP 
SIMULATIONS LISTED IN THE TABLE ARE DRAWN FROM THE SAME 
PARENT DISTRIBUTION. SEE DISCUSSION IN § 5.1. 



h.75Gyr'^ 


p value 


/l.75GBr-^ 


p value 


0.25 


0.01 


0.55 


0.11 


0.30 


0.01 


0.60 


0.11 


0.35 


0.06 


0.65 


0.04 


0.40 


0.10 


0.70 


0.03 


0.45 


0.11 


0.75 


0.01 


0.50 


0.11 







''Mass fraction of stars in the (younger) SSP (age = 1.75 Gyr) in 
the two-SSP simulation, see § 5.1 



formed over an extended period within the cluster itself. 

Three specific scenarios have been proposed for an ex- 
ternal origin of stars with unusual CMDs: (i) the merger 
of two (or more) star clusters, (ii) the merger of a star 
cluster with a giant molecular cloud (hereafter GMC), 
and (in) the entrapment of pre-existing field stars during 
cluster formation. As to the possibility that NGC 1846 
resulted from the merger of two star clusters, the LMC 
is known to host tens of young (10^-10^ yr old) can- 
didate double star clusters (e.g., Bhatia & Hatzidim- 
itriou 1988), many of which have been shown to be 
physically linked objects (e.g., Leon et al. 1999). On 
the other hand, most formation scenarios (Fujimoto & 
Kumai 1997; Theis 2002) suggest that binary clusters 
form within w 10 Myr, much faster than the age range 



seen within NGC 1846. However, it has been suggested 
(Mackey & Broby Nielsen 2007) that star clusters which 
form in groups (star cluster groups or SCGs, see Leon 
et al. 1999) inside GMCs can have ages separated by a 
longer period of time. This idea seems consistent with 
the finding of Efremov & Elmegreen (1998) that the av- 
erage age difference between pairs of star clusters in the 
LMC increases systematically with their (deprojected) 
angular separations, suggesting that the star formation 
process is hierarchical in space and time, i.e., that larger 
star-forming regions form stars over longer periods than 
small regions. However, the very small star-to-star vari- 
ation in [Fe/H] in NGC 1846 allowed by the narrowness 
of its RGB seems difficult to account for in this sce- 
nario: Fig. 1 of Efremov & Elmegreen (1998) shows that 
an average age difference between star cluster pairs of 
~ 200 Myr corresponds to angular separations >0?2, or 
> 150 pc at the distance of the LMC. This size is larger 
than that of any GMC found in the Galactic surveys re- 
viewed in Efremov & Elmegreen (1998). Coupled with 
our finding above that the main sequence turnoff region 
in NGC 1846 does not appear bimodal as much as it ap- 
pears continuous, a binary cluster merger origin for this 
cluster seems unlikely. 

Recent simulations of mergers of star clusters with 
GMCs showed that such mergers can form a second- 
generation of stars that is (at least initially) more cen- 
trally concentrated than the first generation of stars 
(Bekki & Mackey 2009), which is consistent with what 
we find in NGC 1846. However, to reproduce the ob- 
served features in intermediate-age star clusters with 
broad MSTO regions, this scenario seems to require 
rather strongly constrained ranges of GMC parameters 
such as their spatial distribution, mass function, dynam- 
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ics, and chemical composition at the time the mergers 
would have occurred. For example, all currently known 
star clusters in the LMC with broad MSTO regions ex- 
hibit a range in age of ^ 1 — 3 x 10^ yr (cf. Mackey et al. 
2008; Milone et al. 2008b). Larger age differences have so 
far not been found. This would require significant fine- 
tuning of the range of distances between the GMCs and 
the initial star clusters to be merged and/or their relative 
velocities, whereas the current locations of the star clus- 
ters with broad MSTO regions cover a rather wide range 
within the LMC. Another constraint implied by this sce- 
nario is that the GMCs would all have to have the same 
[Fe/H] as the star clusters with which they merged, given 
the very small range in [Fe/H] allowed by the narrowness 
of the RGB in star clusters like NGC 1846. 

Other simulations suggest that during formation, star 
clusters can trap pre-existing field stars, leading to the 
presence of stars with a range of ages (Fellhauer et al. 
2006; Pflamm-Altenburg & Kroupa 2007). This scenario 
also makes predictions which are consistent with observa- 
tions: (i) the radial surface density of the trapped (older) 
stars will be flatter than that of stars formed in the clus- 
ter itself, and (ii) the absence of clear peaks in the distri- 
bution of stars across the MSTO region also seem com- 
patible with this scenario. The biggest point against this 
scenario is that the simulations of Fellhauer et al. (2006) 
indicate that a star cluster with an initial mass of ~ 10 
Mq would only be able to trap up to a few percent of 
its initial mass. This is inconsistent with the 30-40% 
of the stars within NGC 1846 which appear to belong 
to an older subpopulation. The apparent homogeneity 
in [Fe/H] of all stars in NGC 1846 also seems hard to 
explain within the field star trapping scenario. 

We find that a situation where all stars formed in situ 
is in good agreement with several different findings for 
NGC 1846 noted in this paper. In particular, we sug- 
gest a scenario where gas lost by means of low-velocity 
winds of stars with masses lower than those producing 
supernovae (SN) of type accumulates in the central 
regions of the star cluster and forms a second generation 
of stars on timescales similar to the observed age dif- 
ference between the upper and lower end of the MSTO 
in NGC 1846. The stellar agents responsible for such 
self-enrichment may be (a) slow winds of fast-rotating 
O- and B-type stars, which occur at ages of 1-10 Myr 
(Decressin et al. 2007a) or (h) low- velocity winds of 
intermediate-mass (~ 3 7 M©) stars in the AGB phase, 
which have ages of 100 300 Myr (Ventura et al. 200L 
2002; Ventura & D'Antona 2008). Note that the time 
scale for intermediate-mass AGB stars to form in the first 
population is similar to the actually observed age differ- 
ence between the upper and lower end of the MSTO in 
NGC 1846. D'Ercole et al. (2008) recently performed hy- 
drodynamical and A^-body simulations to study the dy- 
namical evolution of star clusters in the context of such 
a self- enrichment scenario. Interestingly, the D'Ercole 
et al. results show that clusters with masses, ages, and 

1° NGC 1846 has an estimated mass of « 2 X 10^ Mq, based 
on an age of 1.7 Gyr, [Fe/H]= —0.50, and either the Maraston 
(2005) or Bruzual & Gharlot (2003) models assuming a Kroupa and 
Chabrier IMF, respectively. This mass is approximately an order 
of magnitude lower than that necessary to retain the chemically 
enriched gas expelled during the earliest stages of the first collapse 
(e.g., Recchi & Danziger 1995; Bastian k, Goodmaji 2006). 



structural parameters similar to those of NGC 1846 are 
expected to show a second-generation population that is 
significantly more centrally concentrated than the first- 
generation stars (see in particular their Figure 18), just as 
we see in NGC 1846. This would also be consistent with 
the observed effect that the stars in the upper part of the 
MSTO arc more centrally concentrated than stars in the 
RGB and AGB (see Fig. 8), since the latter are expected 
to contain stars from both populations. The slow winds 
from either fast-rotating massive stars or intermediate- 
mass AGB stars would also lead to chemical enrichment 
of light elements in the subsequent generation of stars 
due to products from the CNO cycle. If this scenario is 
correct, then one would expect to see significant and cor- 
related variations in the light element abundances (e.g., 
C, N, O, F, Na, Mg) among the stars in NGC 1846, per- 
haps similarly to those found in several Galactic GCs. 
Measuring the chemical composition of RGB stars in 
NGC 1846 and other LMC clusters with wide MSTOs 
(e.g., BerteUi et al. 2003; Milone et al. 2008b) is feasi- 
ble with current spectrographs on 8-lOm-class telescopes 
and should provide additional relevant evidence to help 
decipher the most likely formation scenario. 

6. SUMMARY AND CONCLUSIONS 

We have used deep 5V"/ photometry from HST/ACS 
images to constriict CMDs of the intermediate-age star 
cluster NGC 1846 in the LMC. We have used the ePSF 
fitting technique developed by J. Anderson (Anderson 
& King 2006), which returns high- accuracy photometry 
of cluster stars extending some 5 magnitudes below the 
main sequence turnoff for this cluster. 

The CMD for NGC 1846 shows a number of interest- 
ing features (many of which have been noted previously): 
(i) a very narrow RGB; (ii) an MSTO region that is 
clearly broader than the (fainter) single-star main se- 
quence; (iii) an obvious sequence of unresolved binary 
stars, somewhat brighter and redder than the single star 
main sequence. We have performed a number of new ex- 
periments, including checking the effects of background 
subtraction, completeness corrections, and binary star 
evolution, which confirm that the cluster has a statisti- 
cally significant broad MSTO region. 

We fit isochrones from the Padova, Teramo, and Dart- 
mouth groups in order to determine the best-fit popu- 
lation parameters for NGC 1846. We perform two fits 
with each set of isochrones: One for the upper half of the 
MSTO region and one for the lower half. All three sets 
give a reasonably good fit to the CMD, although there are 
differences between the observations and predictions in 
the shape of the RGB for scaled-solar abundance ratios. 
The overall best fit to the entire CMD is achieved using 
the Dartmouth isochrones. The adopted parameters for 
NGC 1846 are: Distance (m - M)o = 18.45, reddening 
Av = 0.09, ages of 1.7 Gyr (upper haff of MSTO) and 
2.0 Gyr (lower haff), metallicity [Fe/H] = -0.50, and a 
non-solar abundance ratio [a/Fe] = +0.2. 

We use the results from the isochrone fitting to quan- 
tify typical systematic errors of fitted population param- 
eters for intermediate-age star clusters like NGC 1846 
introduced by using any one family of isochrones, the in- 
clusion or exclusion of convective overshooting, and the 
assumption of solar abundance ratios. These systematic 
errors are typically of order 15 - 30% for any given output 
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parameter mentioned above. 

We find for tlic first time that tlic radial distributions of 
stars of different masses/ages vary witli statistical signif- 
icance within NGC 1846, with stars in the upper MSTO 
(i.e., the younger stellar generation) being more centrally 
concentrated than stars in any other area of the CMD, 
including more massive RGB and AGB stars. Since this 
cannot be due to dynamical evolution of a simple stellar 
population, we conclude that the upper and lower MSTO 
regions correspond to intrinsically different populations 
which have undergone different amounts of violent relax- 
ation during their collapse. 

We have tested the role played by binary evolution 
on stars in the MSTO region of the CMD via Monte- 
Carlo simulations of multiple stellar generations. Our 
multi-SSP models include late stellar evolutionary stages 
past the He-corc flash, a realistic treatment of photomet- 
ric uncertainties, and a flat distribution of primary-to- 
secondary stellar mass ratios for binary stars. A quanti- 
tative comparison of the distribution of the stars in the 
MSTO region with those in the simulations that incor- 
porate two SSPs with an age difference as indicated by 
the isochrone fitting shows for the first time that the ob- 
served CMD of NGC 1846 can be better explained by a 
population with a distribution of ages rather than by two 
discrete SSPs as suggested before. 

Finally, we compare predictions of various formation 
scenarios with the observed properties of NGC 1846 de- 
scribed above. Scenarios where some fraction of the 
stars originated from outside the cluster, whether be- 
cause NGC 1846 is the product of a binary cluster merger 
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or a merger of a star cluster with a giant molecular cloud, 
or because the star cluster accreted older field stars dur- 
ing the formation process, do not appear to easily match 
all observed properties. A scenario where a significant 
fraction of the stars form later within the cluster itself 
out of gas enriched by stars of the earliest population(s) 
that feature slow stellar winds appears more consistent 
with the observed stellar radial density distributions and 
the age range implied by the width of the MSTO region. 
Viable sources of the enriched material are thought to in- 
clude fast-rotating massive stars and intermediate-mass 
AGB stars. The latter have ages similar to the observed 
range in ages across the MSTO region. Detailed abun- 
dance ratios from high-resolution spectroscopy of indi- 
vidual cluster stars should be very useful in constraining 
the possible source(s). 
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